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Let /, U and C represent, respectively, the free energy, the internal energy and 

the specific heat of the critical Ising model on the M x N square lattice with periodic 
1 (I , 

boundary conditions, and /oo represents / for fixed M/N and N — ► oo. We find 
that /, U and C can be written as: N(f - /«,) = J2Zi f2i-i/N 2i ~ 1 , U = -\/2 + 
££i u 2i -i/N 2i - 1 and C = 81niV/7r + ££ d/ N \ i- e - N f and u are odd functions 
of N^ 1 . We also find that U2i-i/ 'on-i = 1/V2 and v,2i/c2i = for 1 < i < oo and 
J> . obtain closed form expressions for /, U, and C up to orders 1/iV 5 , 1/iV 5 , and 1/./V 3 , 

respectively, which implies an analytic equation for C5. 
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I. INTRODUCTION 



The Ising model has been used to represent critical phenomena in ferromagnets, binary 
alloys, binary fluids, gas-liquid mixture, etc and is perhaps the most widely studied model 
of critical phenomena 0. For analyzing the simulation or experimental data of finite critical 
systems [0, it is useful to appeal to theories of finite-size corrections [|J and finite-size scaling 
H . Such theories have attracted much attention in recent years @,||J^J|,|] because of the fast 
advance in computers' computing power and algorithms for simulating or analyzing data. 
Theories of finite-size effects and of finite-size scaling in general have been most successful 
in deriving critical and noncritical properties of infinite systems from those of their finite 
or partially finite counterparts. Finite-size corrections and finite-size scaling for the M x N 
square lattice Ising model are of particular interest because the Ising model is very popular 
and such system is usually used to test the efficiency of algorithms for studying critical 
systems [TT5| . In the present paper, we present new analytic results for finite-size effects in 
the Ising model on a large M x N square lattice at the critical point. 

Finite-size scaling is the basis of the powerful phenomenological renormalization group 
method |llJT2| . In the two-dimensional Ising model the finite-size effect on the renormaliza- 
tion transformation has been demonstrated to be rather benign [13|, and the effects due to 



convergence to the fixed point and finite-size are clearly distinguished [l3|. The finite-size 
scaling theory predicts that near the critical point the singular part of the thermodynamic 
quantity of a finite system, say Q s , has the scaling form 



Q a = WY Q (Llt 



(1) 



where L is system linear size, ^ is the correlation length of the bulk system, yg is a 
critical exponent and Yq is the scaling function. The scaling ansatz mentioned above ignores 
the possible logarithmic corrections. In the case of planar Ising model, which displays a 
logarithmic singularities in the specific heat behavior due to a relation between scaling 
exponents in the renormalization group theory pT| , the scaling form ([]]) must then be 
replaced by a more general form p, |IB , nj| , [I7] 



Q s = W Y Q {L/U) + L y Q In L X Q (L/£ C 
which in the case of the specific heat (C) becomes 

C s = Y C (L/U) + ^LX C (L/U)- 



(2) 



(3) 



The results of this paper to be presented below show that the leading term of C s is 8 In L/tt, 
all other finite-size corrections to the specific heat are always integer powers of L" 1 , which 
also imply that the scaling function Xc in Eq. (|3]) is constant and equal to 8/V. Very 



recently, Caselle et al. \TB\ have shown that this result can be predicted by conformal field 
theory under a number of general conjectures. 

The relevance of the finite-size properties to the conformal field theory is another source 
of interest. Discussion of general properties of non- universal corrections to finite size scaling 



and their relation to irrelevant operators in conformal field theory can be found in |]1| . On 
the basis of conformal invariance, the asymptotic finite-size scaling behavior of the critical 



free energy f^ per site and the inverse correlation length ^ N l of a N x oo system is found 
to be || 



lim NiJ = 2ira, (5) 

where /^ is the free energy of the bulk system, c is the conformal anomaly number and 
x is the scaling dimension. The corrections to Eqs. (|]) and (|5|) can be calculated by the 
means of a perturbated conformal field theory |H|,[2l| and can be expressed in terms of the 
universal structure constants (C n i n ) of the operator product expansion [BT|. Quite recently, 
Izmailian and Hu || studied the finite-size correction terms for the free energy and the 
inverse correlation length of critical Ising model on N x oo lattices and obtained a new 
set of the universal amplitude ratios for the coefficients in the free energy and the inverse 
correlation length expansions. It was shown that such results could be understood from a 
perturbated conformal field theory. 

Based on Onsager's solution, explicit calculations of the specific heat finite-size scaling 
behavior have been reported by Ferdinand and Fisher 0] and by Kleban and Akinci ||23||. In 



1969, Ferdinand and Fisher || first studied finite-size corrections for a critical Ising model 
on M x N square lattices with periodic boundary conditions. They gave explicit expressions 
for the critical free energy /, internal energy U, and specific heat C per lattice site for a 
fixed £ = M/N and large N up to orders 1/N 2 , 1/N, and 1/N, respectively: 



f = foo + 



£N 



In (0 2 + 3 + 4 ) -i In (40 2 3 4 ) 



°G^ < 6 > 



^-V2-i^VW4) (7) 



iV 2 + 3 + 4 \N 2 

_ 8 , Ar 8 /, 2 5 / 2 ^ vr\ 
C = - In N+ - In + C S -- 

7T 7T V 7T 4y 

16 Ei 2 t ln0 t / 2 3 04 \ 2 f , 

7T 2 + 03 + 04 V02 + 03 + 04/ 

where 0, = 0j(O, g) (z = 2,3,4) is elliptic theta functions of modulus g = e _?r ^, Ce is the 
Euler constant and f^ is the free energy in the thermodynamic limit M, N — > oo. 

In 1983, Kleban and Akinci [^J gave a very accurate and relatively simple approximate 
closed form expression for leading specific-heat correction term which results from retaining 
only the two largest eigenvalues of the transfer matrix. This approximation is already good 
at £ = 1 and becomes exponentially better with increasing £, and they interpreted their 
results in terms of domain-wall energies. 

In this paper we study the same system as [H] and find that /, U and C can be written as: 
N(f-foo) = J2Zi f2i-i/N 2i ~\ U = -V2+E,=i ^-i/iV 24 " 1 and C = 8 InN/n+EZo d/N\ 
i.e. Nf and U are odd functions of A r_1 . We also find that u^i-i/ C2i-i = l/v2 and 



U21/C21 = for 1 < % < oo and obtain analytic equations for /, U, and C up to orders 1/N 5 , 
1/N 5 , and 1/iV 3 , respectively, which implies an analytic equation for C5. 

We have also shown that Kleban and Akinci approximation is in excellent agreement 
with our exact results for the leading correction terms of the free energy (f\), the internal 
energy (u\) and the specific heat (co,Ci) at £ > 1. For the next correction terms the error 
introduced by the two-eigenvalues approximation is maximum at £ = 1 (M = N). With 
increasing £ the exact and approximate values approach exponentially and approximation 
becomes already good at £ = 1.65 for the correction terms f 3, 113,02, C3 and at £ = 1.85 for 
the correction terms $5,1x5. 

This paper is organized as follows. In Sec. II, we write the free energy /, the internal 
energy U , and the specific heat C of the Ising model in terms of P±, Pi, P%, Pi, Qi, Q2, and 
Qs defined in this section. In Sec. Ill, we present asymptotic expansions for /, U, and C. 
In Sec. IV, we discuss some problems for further studies. Some mathematical details used 
in the derivation of equations in Sec. Ill are given in Appendix A. 

II. ISING MODEL 



Consider an Ising ferromagnet on an M x N lattice with periodic boundary conditions 
[i.e. a torus). The Hamiltonian of the system is 



^ = -^E 



S;S 



l-JJI 



(9) 



<ij> 



where /3 = (/c^T) -1 , the Ising spins Sj = ±1 are located at the sites of the lattice and 
the summation goes over all nearest-neighbor pairs of the lattice. The partition function 
Zmn(T) of a finite M x N square Ising lattice wrapped on a torus can be written as 

1 



Z MN {T) = \{2smh2J)^ MN Y,Z l , 
where the partial partition functions Zj are defined by 



N-l 



M 



Zi = Yl 2 cosh — 7 2r +i = Pi exp 

r=0 ^ 

"- 1 M 

Z 2 = [[ 2 sinh— 7 2r +i = P 2 exp 

r=0 ^ 



~7T 2^ ^2r+l 

~7T 2^ 72r+l 
Z r=0 



with 



^r 1 , M 
Z 3 = [[ 2 cosh— 7 2r = 

r=0 ^ 


P 3 


exp 


\M N ^ 1 ' 

■9- Z^72r 
. Z r=0 


(l + e-^-w) 


N-l M 

z 4 = n 2 sinh v 72r = 

r=0 Z 


Pi 


exp 


"M^ 1 1 


(1 - e- M7 °) 


N-l 

P 1= n(l + e' M72 " +1 ). 

r=0 






N-l 

P 2 = J] (1 - e - M ^ +1 

r=0 


JV-1 

p 3 = n( i + e_M72r )> 

r=l 






JV-1 
r=l 


( 1 _ e -M 7 2.) 



(10) 

(11) 

(12) 
(13) 
(14) 



(15) 



and 7 r is implicitly given by 



cosh 1J r-a 

cosli7 r = , - - - — cos — -. 
' sinh2J N 



(16) 



At the critical point J c of the sq lattice Ising model, where J c = |ln(l + v2), one then 
obtains 



7^ cr) = 2ip sq ( — j , with ip sq (x) = \n(sinx + \fl + sm 2 x). 



(17) 



The free energy, the internal energy per spin and the specific heat per spin can be obtained 
from the partition function Z MN 



1 1 1 1 / 4 \ 

/ = — — In Z MN = - In 2 sinh2J + — — In- V Z t ), 
J MN 2 MN 2 \{r{ % ) ' 

> MN = -coth2J- J^ (j>;) / (j>) > 
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1 d 



MN dJ 
1 d 2 



MNdJ 2 

2 



\nZ MN 
1 



(18) 
(19) 
(20) 



sinlT 2 J MN 



+ SF E^/E^ 



\i=l 



vi=l 



E^ / £*< 



vi=l 



u=l 



where the primes denote differentiation with respect to J. At the critical point (T = T c 
the partial partition functions Zi and their first and second derivatives are given by 



R 



Z l = P l e A , Z 2 = P 2 e A , Z 3 = 2P 3 e", Z 4 = 0; 

Z[ = 0, Z' 2 = 0, Z 3 = 0, Z; = 4MF 4 e B ; 

yll yll yll 

Ql^li nr-nT = Q2Z2, TT1TF = Q3Z3, Z 4 = —y2Z 4 , 



MN 



where 



MN 



A _M N ' (cr) 

n — 9 Z_^ 72r+l) 



MAT 



M^ 1 (cr) 

Z r=0 



(21) 
(22) 

(23) 



(24) 



^ J- v-^ "( CJ *H l, M 72r+1 

Gi = 2N E ^r+i tanh , 



JV-1 



(cr) 



/-) 1 V^ "( cr ) +1, M 72r+1 
^2 = 77^7 E 72r+l COth -^-. 



2N 



r=0 



JV-1 



"(«-) + „„l, M T: 



^3 = 4e + ^E72rtanh 



2r 



(25) 
(26) 

(27) 



and Pj are given by Eq. (|15D with j r = 7^ and 7^ cr ^ denote the second derivative of 7,, 
with respect to J at the critical point J = J c . Then the exact expression for the free energy, 



the internal energy and the specific heat of a finite Ising model at critical point (T = T c ) 
can be written as 

/ = lh,2 + J-A + JL h ,fl + ft + 2ft« 1 >(i>-A) 

^ 2 MAT M^V 2 ' v ; 

7/ = _*/2 — — — (29) 

N2P 3 + (P 1 + P 2 )exp(A-By K } 



c = V2U -n 



4P 4 



2 



2P 3 + (Pi + P 2 ) exp (A - B) 



2Q3P3 + (Q1P1 + Q2P2) exp (A - B) 

2P 3 + (P 1 + P 2 )ex P (A-B) • [ ' 



III. ASYMPTOTIC EXPANSIONS 

We consider only sequences of lattices in which £ = M/N remains positive and finite as 
the thermodynamic limit M, N —>■ 00 is approached. Using Taylor's theorem we find that 
Mj r is even function of 1/N at the critical point 



x <H vr 3 £ r 3 vr 5 £ r 5 

iV^ ~ ^ r ~ UJP + ~96~ TV 4 



M 7 -) = y ^- = ^r - ^— + ^— + .... (31) 

lr Z_j AT2i s 1 O /\r2 Ofi AH K ' 

8=0 



Using Euler-Maclaurin summation formula |24j] we can expand A and -B up to arbitrary 
order 

2G_ AT tt£ 7tt 3 £ 1 31vr 5 £ 1 10033vr 7 £ 1 , , 

— MN A - -I - I - I - I- (32) 

it 12 1440 N 2 24192 N 4 9676800 JV 6 '"' l ; 

2G„ rA7 tt£ tt 3 £ 1 7r 5 «e 1 79tt 7 £ 1 

= MN - - - - (33) 

■k 6 180 N 2 756 N 4 75600 A^ 6 '"' y ' 

where B 2 i are the Bernoulli numbers and G = 0.915965... is Catalan's constant. 

Let us now evaluate the products Pi for i = 1,2,3,4. It is easy to see from Eqs. fllBD 
and ( |31| ) that Pi contains only even power of 1/N 

"-(O.Oli + E^I (. = 1,2,3,4), (34) 

«=*(o,e) (i + £ + *! + ...), * =*«u)(i + £ +£+... 

P 3 = ,s(0,« 1 4- || 4- f§ 4- ...) , P 4 = , 4 (0,fl 1 4- |i 4- £ 4- ... 



with 



7^(0,0 = ^, 7T a (O,0 = £, 7T 3 (0,0 



_0j_. 

9^9 A 



vr 4 (0,0 = ^o 5 



e * 



20 o 3 
9 2 9 3 9 4 



where 0« = 0j(O,g) is elliptic theta functions of modulus q = e~^. The explicit expressions 
for the coefficients pn and pi2 for i — 1, 2, 3, 4 are given in Appendix A. 

One readily sees from Eqs. (pip, (|29| ) and fl3~2|) - fl3"4|) that the finite-size estimates of the 
free energy (Nf) and the internal energy (U) must be odd functions of A" -1 . 



A (,/ — Jooj = 2^ jy2i-l ' 



i=l 



^-v/2 + ^r 



«=i 



«2i~l 

]\T2i-l ' 



(35) 
(36) 



Substituting Eqs. (||) - ([34]), (TAT4D and (TAT5D in Eqs. © and (g9|) we can write the 
expansions of the free energy (Nf) and the internal energy (U) at the critical point (T = T c ) 
up to 1/A 5 order. The final result is 



with 



and 



with 



s c 



-1 r 



a(/ - /oo) = y 



In (0 2 + 3 + 04) 



-In (40 2 0304) 



7T 3 80 2 0304 [0 3 3 (0 2 3 + 01) 



0f0l] 



7(0 2 9 + 3 9 



04 9 ) 



A 3 



+ 



1440(02 + 03 + 04 

+ 



1 £A 8 hi + (#-#) /52 , ^ / J_ 



A 5 1897T 2 2 + 3 + 04 



VA 7 



/ 52 = 2 (-32 + 48A; 2 - 78A; 4 + 31A; 6 ) + 3 (31 - 78k 2 + 48k 4 - 32k 6 ) 

+ 4 (31 - 15k 2 - 15A; 4 + 31k 6 ) 
f sl = 2 (-32 + 80k 2 - 38k 4 + 21k 6 ) + 3 (31 - 88k 2 + 88A; 4 ) 

+ 4 (31 - 67k 2 + 25k 4 - 21k 6 ) 

21 3 4 A; 4 (1 + k' 2 ) 2 + 2 3 A;' 4 (1 + k 2 ) 2 + 2 4 (A;' 2 - k 2 ) 2 



+ 



U 



V2- 



020304 



02 + 03 + 04 

2 020304 (0 2 9 + 0f + 0f) vr 3 e 



+ 



+ 



A 2 + 3 + 4 ' A 3 

020304 e 2 ^ 8 



A 5 (0 2 + 3 + 4 ) 2 3vr 2 



u 51 + 



h + 03 + 04) 2 

'A' A' 



AT' A 



U, 



r ,2 



96 



1 

aV 



^ 



^52 

t/51 



3 (A;' 2 - k 2 ) - 9 2 k 4 (l + k' 2 ) + 9 4 k'\l + A; 2 



0. 



23 - 64A; 2 + 64A; 4 



9,k 4 



16 + 8A; 2 - k 4 



+ 9Jl-k 



2\2 



+ 



24 24 

3 4 A: 4 (1 + A;' 2 ) 2 + 2 3 A: /4 (1 + A: 2 ) 2 + 2 4 (A;' 2 - A; 2 ) 2 
12(02 + 03 + 4 ) 

7 



23 - 6A; 2 - A; 4 
24 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



where /oo = —0.5 In 2 — 2G/n and 82, #3, #4 are elliptic function 



*=,/«, 8, = ^, ft = ^/™ (43) 

with K(k) and £/(fc) the elliptic integrals of the first and second kind, respectively. For 
simplicity we denote K = K(k), K' = K'(k), E = E{k) and E' = E'(k). 

It is interesting to compare our results for the free energy and the internal energy with 
Kleban and Akinci two eigenvalues approximation. Keeping only two largest eigenvalues Ao 
and Ai of the transfer matrix, the partition function of the Ising model can be written as 

Zmn = A* 7 + Af (44) 

with 

/-, N-l \ 

A = (2 sinh 2J) N / 2 exp f - £ 72r+1 j , (45) 

Ai = (2 sinh 2J) N ' 2 exp Q £ l2 \ , (46) 

where 7^ is implicitly given by Eq. (]T6|) . 

To write the critical free energy / and critical internal energy U in the form of Eqs. ([55]) 
and (|36|), we must evaluate Eqs. fl45|) and (|46| ) asymptotically. These sums can be handled 
by using the Euler-Maclaurin summation formula |24]|. After a straightforward calculation, 
we have obtained 

x(<w) Aapp) Aapp) 

w-/oo) = v + V + V + -' (47) 

„(«PP) „,(app) ..(app) 

l, =-^ + V + V + V + - (48) 

with 

/&•*> = -Jl - 1 In (l + e -*/*) (49) 

/3(aPP) = 2^0 [1 " 15tanhK/8)] (50) 

^ = 48^4 [1 " 63tanh (^/8)] - ^sech 2 K/8) (51) 

y<fiPP) = _i + tanh (tt£/8) (52) 

4 W) = ^sech 2 «/8) (53) 

4^ = ^|sech 2 K/8) - |^sech 2 K/8) tanh K/8) (54) 

The expressions of the coefficients given by Eqs. (f!9|) - ( |54"D are much simpler than their 
exact counterparts given by Eqs. ([37]) - (|42|) . Nevertheless, one can see from Figs. 1 and 2 

8 



that two-eigenvalues approximation proposed by Kleban and Akinci is already good at £ = 1 
for the leading corrections terms in the free energy (fi) and the internal energy (ux) and 
becomes exponentially better with increasing £. The error introduces by the two-eigenvalues 
approximation is maximum at £ = 1 (M = N) . With increasing £ the exact and approximate 
values approach exponentially and approximation becomes already good at £ = 1.65 for the 
correction terms f'3,11,3 and at £ = 1.85 for the correction terms /s,U5. We consider the case 
£ > 1 only. By symmetry, the same results hold for £' = l/£ < 1. 

To calculate the specific heat we must also evaluate asymptotically the sums appearing 
in the expression fl3Tf ) for C, namely Q±,Q2 and Q 3 . Since the analysis follows the same 
general lines as in the cases of the free energy and the internal energy, we will not present the 
details of calculations and we quote here only results, namely, at the critical point T = T c 
the asymptotic expansion of the sums Q\, Q2 and Q3 can be written as 



Q, = -lniV + ]r 



7T 



3=0 






for 



1,2,3) 



where q i0 and q a (for i = 1, 2, 3) are given by 

8 2 5 / 2 
q w = -(C E + ln 2 In 3 ) 

71 71 

8 2 5 / 2 
g 2 o = -{C E + ln 2 ln0 4 ) 

71 71 



2 5/2 



<?30 = -(C B + ln 

7T 71 

9n = 



21n# 2 ) 



^21 
931 



9tt 2 
8^^ 
~ 9tt 2 
8K 4 £ 
9tt 2 



— 2 )(f-r 



1 - 3fc 2 + (1 + k 2 



' E' E' 



2-se + (2-e) (--- 



K> K 



(55) 



(56) 
(57) 
(58) 
(59) 

(60) 

(61) 



It is easy to see from Eqs. fl29]), (|30|), (0) and (|55| ) that the asymptotic expansion of the 
specific heat, can be written as 



C = -\nN + y Cl /N\ 

71 n 



(62) 



Except for the leading term, all other corrections in the asymptotic expansion of the spe- 
cific heat are proportional to 1/N 1 , without multiplicative logarithms. This result imply 
immediately that scaling function Xc in Eq. @ is constant and equal to 8/tt. 

It is also clear that the contribution to odd (iV~ 2 * -1 ) order in the specific heat expansion 
give only first term in right-hand side of the Eq. (|30|). Thus, we can obtained immediately 
that the ratio M2«+i/c2j+i of subdominant (iV -2 * -1 ) finite-size corrections term in the internal 
energy and the specific heat expansions are constant, namely, 



U2i+l/c 2 i+l = l/\/2 



(63) 



as well that w 2 j/c 2 j = for 1 < i < oo. 

Let us now evaluate the first few terms in the specific heat expansion. Substituting Eqs. 
(0) - (H), (§§, d|) - (|ip, (|M§ and ( pQ5| ) in Eq. Q3IJ) we have finally obtained the 
expansion of the specific heat (C) at the critical point (T = T c ) 



c= s lnJV+ s( ln ^ + c *)_ 4 £/ — 



2 



16E-=2^1n^ n ^ 2 3 04 1 c 2 



2^^^ + -^ (64) 



7T 2 + fl 3 + #4 02 + 03 + 04 N N 2 

N*24V2 (02 + 03 + 0a) 2 ViVV + \N* 



with 



71^ 010WM + 01 + 01) ^ vr 2 £ 0I0\{20 2 -0 3 - 4 ) 

C 2 — ~ T^ 77 7~?) , n \q r 



12 (# 2 + #3 + # 4 ) 3 9 # 2 + # 3 + #4 

7T 2 e 2 #3#4 



^ + 0t)0l ^ ^ - (^ + *M ^ £ + (^ - *M ^ £ 

#4 #3 #4 



6 (fl 2 + # 3 + #4) 2 
* gj + flg + W 2 4 + <%) ~ Z0 2 03{01 + 01) A. 9 , fl2 p \ rfi ,^ 

Equation ( |63| ) imply that the amplitude of the term 0(1/N 5 ) in Eq. (|64]), i.e. C5, is v^s 
where W5 is the amplitude of the N~ 5 correction terms in the internal energy expansion Eq. 

In two-eigenvalues approximation the specific heat can be written as 

o „( a PP) J a PP) Jo-PP) 

c = ^i„iv + 4°» ) + V + V + V + - (66) 

with 

c { app) = - fin — + C E - -\ + £ sech 2 (7r£/8) + — [-1 + tanh (tt£/8)] (67) 

7T \ 7T 4 / 7T 



C 



\/2[-lH-tanh(7T^/8)] (68) 



(app) 
1 

4 aw) = -| + I hi + tanh fa/8)] + — ^sech 2 K/8) 

7T 3 ^ 



96 
3, 



sech 2 «/8)tanh(7r^/8) (69) 



4 aw) = ^|sech 2 «/8) (70) 

We plot the aspect-ratio (£) dependence of the finite-size specific heat correction terms Cq 
and c 2 in Fig. 3. The exact and approximate values approach exponentially as £ increases. 
Note, that the ratios of correction terms U\jc\ and M3/C3 are constant and given by Eq. d63|). 
In Fig. 4 we plot the aspect-ratio dependence of the error introduced by two-eigenvalue 
approximation for the correction terms in the free energy, internal energy and specific heat 
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asymptotic expansions. The deviation of the two-eigenvalues approximation from exact 
result is about one percentage at £ = 1 for the leading correction terms f\, U\, c , at £ = 1.65 
for the second correction terms / 3 ,m 3 , c 2 , at £ = 1.85 for the third correction terms f^,u 5 
and diminishes very rapidly as £ increases. 

It is of interest to compare this finding with other results. Equations ([J7|), (|40|), and ( |64"D 
are consistent with Ferdinand and Fisher's similar expansions || up to orders 1/N 2 , 1/N, 
and 1/N, respectively. Others terms in our equations, except the term of 0(1/N 3 ) for U 0, 
are new. For £ = 1, we have u 3 = 0.206683145 . . . and u 5 = 0.730182312347. . . which are 
quite consistent with numerical data u 3 = 0.206683133 and u 5 = 0.73018231235 obtained 
by Salas and Sokal M. 



IV. DISCUSSION 



The results of this paper inspire several problems for further studies: (i) can one obtain 
an exact asymptotic expansion for the thermodynamic functions up to arbitrary order, as it 
can be done for the Ising model onAfxoo square, honeycomb, and plane triangular lattices 
||. (ii) It is of interest to know whether the amplitude ratio of Eq. ( J63|) can be extended to 
honeycomb and plane triangular lattices, i. e. whether the ratio is universal, (iii) If, so, how 
do such amplitudes behave in other models, for example in the three-state Potts model? 

Note added: After the completion of this paper, we learned that similar results have been 
independently obtained by Salas 121 
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APPENDIX A: 



Let us now evaluate the coefficients pn and p i2 for i = 1,2,3,4. After little algebra, 
following the general lines of the Ferdinand and Fisher paper [R| , we can obtain the following 
expression for the coefficients pn and p-a 
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P31 






1 %^ ,v, cosh7r 6 , 3 v-/ .NjCOshvrO 
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4 j=1 sinh 2 tt£j 2 j=1 sinh 4 ir£j 

\)o 3 ~ f-lV 



(Al) 



E- 



3 ~ > 

+ oE- 



i=1 sinh vr£j 2 j=1 sinh vr£j 
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j=1 sinh 2 7r£j 2 j=1 sinh 4 n£j 



Ps2 = -Pl 1 + 81 
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47T 6 e : 

81 

47T 6 f 

81 

4tT 6 £ 3 
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(A2) 



with 



1 ~ . . / 1 60 120 ' 

128 ~[ \sinh 2 xj sinh 2 xj sinh 6 xj / 

1 ^ / 1 60 120 N 



^ 3 (x 



3=1 \ 

1 ~ . / 1 60 120 

= > cosh xj 5 1 2 1 « — 

128 • \ ysinh xj sinh xj sinh xj 

I5W ^ + 15 15 ^ 



sinh 2 xj sinh 4 xj sinh 6 xj 



(A3) 



J_ V ( 2 15 15 

16 • J \sinh xj sinh xj sinh xj 

Let us now introduce the following notation: 

00 -1 

sjx) = y 

~[ sinh™ xj 
Then the coefficients pa and pj2 can be rewritten in 



for n = 2,4, 6. 

the more symmetrical way 



(A4) 
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|i2(e/2) - ^i?(0 + 2i2(2£) 



P11 = 3^ 

ft! = ivr 3 e [i?(0 - |i2(e/2) 
P3i = ^ 3 e[ J R(0-2i?(20] 
P4i = ~\^R(0 



(A5) 



and 



* 2 (0 = -*(0 + ^(^/ 2 ) 

* 3 (0 = 2*(20 - *(0 



(A6) 
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with 

R(x) = S 2 (x) + ^S A (x) (A7) 

*(s) = 4 (2S 2 (x) + 155 4 (x) + 155- 6 (x)) (A8) 

lb 

Thus we have shown that the coefficients pn and pi 2 can be expressed in terms of the only 
object, namely S n (x) for n = 2,4,6. The S 2 (x) is given by (see [G5) p. 721) 

1 2(2 -A; 2 

- -| — 

6 3ir 2 n u 



S *( x ) = *+ V _, J K 2 (k) - —K{k)E{k) (A9) 



where x = nK'(k)/K(k) with K(k) and E(k) the elliptic integrals of the first and second 
kind, respectively. The S^x) and S 6 (x) are calculated to be 



90 9tt 2 R2{k) + 3^ 



^0*0 = -7^ - -^J^K^k) + _K(/c)£(/c) (A10) 



+ 8(1 ' fc2 4 +fc4) ^) 
45vr 4 V ; 

_ 8(l-fc 2 + fc 4 ) 4 _ 32(2 - 3k 2 - 3fc 4 + 2fc 6 ) 6 

45vr 4 l ; 945tt 6 l ; 



Thus we are now in position to evaluate R{x) and \l/(a;) given by Eqs. ( |A7| ) and (AS) 
respectively. The result is 



1 4(1 - k 2 + K 4 ^ 
60 + 15tt 4 



*(*) = -- + v 1C : y ^ 4 (fc) (Ai2) 



*( x ) = J_ _ 2(2 3fc 2 3fc 4 + 2fc 6 ) K6 
v ; 1008 63vr 6 v ; v ; 

The expressions for R(2x),R(x/2) and *£?(2x), ^?(x/2) can be written as function of the 
modulus k by using properties of the elliptic functions. 

Thus for the coefficients pn (for i — 1, 2, 3, 4) we have finally obtained 

7^Z (j + atf-MX 

mi 1440 90tt V ; ' 

^ 21 1440 90tt K h K J 

^ (-8 + 8^+7^ 
^ 31 180 90vr V ; ' 

vr 3 e 4(1 - fc 2 + K 4 )e ^ 
^ 41 180 45tt v ; 

After little algebra the expressions for pi 2 (for i — 1, 2, 3, 4) can be written as 
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where for simplicity we denote K = K(k), K' = K'(k), E = -E(A;) and £" = E'(k). 
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FIGURE CAPTIONS 



FIG. 1. Finite-size free energy correction terms (a) j\, (b) f 3 and (c) / 5 as functions 
of the aspect ratio £, which are defined by Eqs. (37), (38), (39), and (43). Solid curves: 
exact values; dashed curves: two-eigenvalue approximations of Eqs. (49)-(51). The exact 
and approximate values approach exponentially as £ increases. 

FIG. 2. Finite-size internal energy correction terms (a) u±, (b) u^ and (c) u§ as functions 
of the aspect ratio £, which are defined by Eqs. (40)-(43). Solid curves: exact values; dashed 
curves: two-eigenvalue approximations of Eqs. (52)-(54). The exact and approximate values 
approach exponentially as £ increases. 

FIG. 3. Finite-size specific heat correction terms (a) cq and (b) C2 as functions of the 
aspect ratio £, which are defined by Eqs. (64) and (65). Solid curves: exact values; dashed 
curves: two-eigenvalue approximations of Eqs. (67) and (69). The exact and approximate 
values approach exponentially as £ increases. 

FIG. 4. The error introduced by two-eigenvalue approximation for the correction terms 
in the free energy, internal energy and specific heat asymptotic expansions, (a) The error 
for the leading correction terms j\, U\ and Cq as functions of £. (b) The error for the second 
correction terms fa, u 3 and C2 as functions of £. (c) The error for the third correction 
terms /s and u§ as functions of £. The vertical axes represents the error defined by Error = 
[a exact — a a PP rox ^ /a exact , where a stand for the correction terms in the free energy, internal 
energy and specific heat asymptotic expansions. Solid curves: free energy; dashed curves: 
internal energy; dot-dashed curves: specific heat. The deviation of the two-eigenvalues 
approximation from exact result is about one percentage at £ = 1 for the leading correction 
terms, at £ = 1.65 for the second correction terms, at £ = 1.85 for the third correction terms 
and diminishes very rapidly as £ increases. 
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